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Abstract. A survey of helminths associated with terrestrial slugs focusing on the invasive Arion vulgaris and the native A. 
ater was conducted on populations from France, Germany, Netherlands, Norway and Poland. In total, 648 terrestrial slugs 
were collected from 18 sample sites, and identified by means of morphological examination, dissection of genitalia and 
molecular analysis using mitochondrial DNA. In addition to A. vulgaris and A. ater, also A. vulgaris/A. rufus hybrids and 
A. ater/A. rufus hybrids were collected. Helminth species were identified based on morphological features and sequencing 
of the 18S and ITS rDNA regions. The parasites included four nematode species: A/loionema appendiculatum, Angiosto- 
ma sp., Phasmarhabditis hermaphrodita, Entomelas sp., two trematode species: Brachylaima mesostoma, Eurytrema sp., 
and one cestode (tapeworm) species: Skrjabinia sp. Alloionema appendiculatum was the most common helminth in the 
investigated slug populations. Furthermore, we found higher prevalence of trematodes in the invasive A. vulgaris compa- 
red with the native A. ater, while differences in the prevalence for nematodes were not as clear. 
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INTRODUCTION 


Parasitism plays an important role in the ecology and 
evolution of terrestrial gastropods (slugs and snails) influ- 
encing the evolution of sexual reproduction, life-history 
traits as well as host resistance leading to host-parasite 
co-evolution. Parasites have a direct impact on life cycles 
of their hosts, and the effects of parasites can be modu- 
lated by environmental factors. Several studies have been 
conducted on the impact of climate change on crop pests 
in relation to natural enemies such as insect predators, 
parasitoids, pathogenic microorganisms, and helminth 
parasites (Gerard et al. 2013; Wilson et al. 2015). The 
presence of helminths in slugs can be influenced by the 
size, age and the spatial 1solation of the host population 
and by habitat characteristics (Baur & Baur 2005; Gerard 
et al. 2013; Wilson et al., 2012; Wilson et al. 2015). 
Different groups of helminths can be associated with 
terrestrial gastropods, but slugs as hosts have been given 
most attention (Ross et al. 2010a; Ross et al. 2010b; Ross 
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et al. 2016). Currently more than 25,000 species of nem- 
atodes have been described, of which around 3,500 are 
parasitic nematodes of invertebrates (Laznik et al. 2010). 
Nematodes parasitise both slugs and snails, however slugs 
are parasitised more frequently and by a greater diversity 
of nematodes than snails (Mengert 1953). This 1s because 
slugs usually inhabit soil, thus increasing their exposure 
to nematodes (Morand et al. 2004; Ross et al. 2010a). 
Currently, representatives of eight nematode families are 
known to be associated with terrestrial slugs: Agfidae, AI- 
aninematidae, Alloionematidae, Angiostomatidae, Cos- 
mocercidae, Diplogasteridae, Mermithidae, and Rhab- 
ditidae (particularly the genus Phasmarhabditis) (Ross 
et al. 2017). These families are known to form a number 
of different relationships with slugs, including parasitic 
(specialist or generalist), phoretic and necromenic asso- 
ciations (Ross et al. 2010b; Ross et al. 2017). Moreover, 
trematodes, mites, sporozoa, ciliates, and cestodes have 
also been described to interact with slugs (Stephenson & 
Knutson 1966; Baur & Baur 2005). 
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Terrestrial gastropods are considered to be one of the 
most successful and diverse animal groups in terrestrial 
ecosystems (Barker 2001). Many of them have become 
Invasive species in the context of expanding their range 
of distribution and generating economic damages (e.g., 
Lissachatina fulica Bowdich, 1822 and Deroceras retic- 
ulatum Müller, 1774) (Hammond & Byers 2002; Ross 
et al. 2010a). One of the 100 most invasive species in Eu- 
rope is Arion vulgaris Moquin-Tandon, 1855, commonly 
known as the Spanish slug. It has probably been unin- 
tentionally introduced into new habitats via plant matter, 
packaging, and waste materials (Kozlowski 2007; Hat- 
teland et al. 2013; Zajac et al. 2017). It is a major defoli- 
ator of plants and causes severe damage in orchards and 
gardens as well as in crops (Gren et al. 2009). This slug 
has been assumed to originate from the Iberian Peninsula 
and spread into Central Europe in the 1950s (Frank et al. 
2002), although recent studies suggest a more northern 
origin, possibly in France (Hatteland et al. 2015; Zemano- 
va et al. 2016). Monitoring the spread of A. vulgaris is 
difficult because the pest 1s morphologically similar to 
the other closely related, large arionids (A. ater Linnaeus, 
1758; A. rufus Linnaeus, 1758; A. magnus Torres Min- 
guez, 1923; A. lusitanicus Mabille, 1868; A. flagellus 
Collinge, 1893) that occur in Europe. Arion ater, A. rufus 
and A. vulgaris can hybridise with each other (Roth et al. 
2012; Dreijers et al. 2013) and introgression has readi- 
ly been shown, especially between A. ater and A. rufus 
(Hatteland et al. 2015; Zemanova et al. 2017). Moreover, 
A. vulgaris may outcompete native slug species because 
of its large size and high population densities (Frank 
2003). Native slugs like A. ater and A. rufus have been 
observed to decline and/or disappear in areas colonized 
by invasive slugs such as A. vulgaris in continental Eu- 
rope (Falkner 1990) and Scandinavia (B. A. Hatteland 
pers. obs.) and a similar pattern seems to occur where 
A. flagellus has been introduced in Britain (Davies 1987). 

There are numerous hypotheses regarding what makes 
species invasive. Release from natural enemies 1s regard- 
ed to be an important factor supporting invasiveness by 
many organisms (Torchin et al. 2003). The enemy release 
hypothesis (ERH) states that the lack of natural enemies 
in an invader's introduced range influences its abundance 
or impact (e.g., estimated using individual size, popula- 
tion abundance, or propensity to displace native species) 
(Torchin et al. 2003; Colautti et al. 2004). Torchin et al. 
(2003) studied parasite burdens in 26 species of invasive 
animals and found that most of them had fewer parasites 
in their introduced areas compared with their home rang- 
es. Indeed, it is less likely that hosts will spread parasites 
into their introduced range since introduced populations 
often originate from relatively small subsets of native 
populations (and sometimes from uninfected life-histo- 
ry stages) (Torchin et al. 2003). However, recent studies 
have shown that A. vulgaris hosts a range of parasites in 
introduced areas as well as relatively high parasitic loads 
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compared with other native slug species (Ross et al. 
2010a; Ross et al. 2016). 

This study describes results on the diversity and dis- 
tribution of helminths associated with A. vulgaris and 
A. ater in Europe, Le France, Germany, Netherlands, 
Norway and Poland. 


MATERIAL AND METHODS 
Collection and identification of slugs and helminths 


Slugs were collected from 18 sites in Europe (France, 
Germany, Netherlands, Norway, Poland) in 2015 and 
2016 from late August to October (Table 1). Sites were 
selected based on information from local growers and 
gardeners as well as advisory services in the region re- 
garding the presence of slugs. At each site more than 10 
slugs were collected and sent to the Institute of Environ- 
mental Sciences (Jagiellonian University, Krakow, Po- 
land) and Norwegian Institute of Bioeconomy Research 
(NIBIO, As, Norway). Slugs were identified by means 
of morphological examination guided by von Proschwitz 
(2009), dissection of their genitalia and molecular anal- 
ysis using a fragment of mitochondrial cytochrome c 
oxidase subunit I (COI, mtDNA). The main features of 
the genitalia in A. vulgaris are the small atrium, almost 
symmetrical, one-partite, bursa copulatrix oval and a 
free oviduct with a short, thin posterior end and a thick, 
rapidly expanding anterior end with a large, asymmetric 
ligula inside (Wiktor 2004). The main characteristics of 
A. ater are the atrium and vagina considerably narrower 
than spermatheca, oviduct narrow and the spermatheca 
spherical (Welter-Schultes 2012). During the dissection 
of slugs, a piece of tissue was taken for DNA extraction 
and preserved in 96% ethanol at -80°C. 

Slugs were checked for potential helminths, which 
were identified using a combination of morphological 
and molecular techniques. The helminths were first clas- 
sified morphologically as nematodes or trematodes un- 
der a stereomicroscope Olympus SZX10. For purposes 
of molecular identification, all helminth samples (e, 
adults, juveniles and cysts) were transferred to Eppen- 
dorf tubes containing 70% ethanol. Each helminth was 
introduced to a separate Eppendorf tube. 


DNA isolation, amplification and sequencing 


DNA extraction of slugs was performed using a commer- 
cial DNA extraction kit (NucleoSpin® Tissue, Mache- 
rey-Nagel, Duren, Germany), which uses a proteinase 
K to digest proteins within cell membranes and columns 
with silica membranes and buffers for cleaning extracted 
DNA. PCR reactions were performed to obtain a frag- 
ment of mitochondrial DNA with LCO1490/HC02198 
primers (Folmer et al. 1994). A PCR reaction was per- 
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formed in a reaction mixture of 20 ul per each sample and 
consisted of 3 ul of template DNA, 0.6 ul of each primer, 
2 ul of 10x buffer, 13 ul of ddH2O, 0.6 ul of 20 mM 
dNTP and 0.2 ul of DreamTaq!M DNA Polymerase 
(Thermo Fisher Scientific Inc., MA, USA). PCR condi- 
tions included 5 min initial denaturation at 94°C and then 
1 min denaturation at 94°C, 1 min 30 s annealing at 45°C 
and 1 min 30 s elongation at 72°C for 5 cycles and then 
1 min denaturation at 94°C, 1 min 30 s annealing at 50°C 
and 1 min elongation at 72°C for 35 cycles followed by 
a final elongation step for 5 min at 72°C. A volume of 5 
ul sample of PCR product was run on a 1.0% agarose gel 
for 30 min at 100 V to check DNA quality. PCR products 
were cleaned up by commercial kit (NucleoSpin® Gel 
and PCR Clean-up, Macherey-Nagel, Duren, Germany). 
The sequencing reaction was performed in a reaction 
mixture of 10 ul per each sample and consisted of 2 ul of 
template DNA, 5.85 ul of ddH2O, 0.15 ul of primer, 1 ul 
of 5x buffer and 1 ul of Terminator (BrightDye® Termi- 
nator Cycle Sequencing Kit, MCLAB, South San Fran- 
cisco, USA). Sequencing products were cleaned by using 
a kit to remove terminators after sequencing reactions 
(ExTerminator, A&A Biotechnology, Gdynia, Poland). 
The sequencing reactions were performed in Molecular 
Ecology Lab, Institute of Environmental Sciences (Jagi- 
ellonian University, Kraków, Poland). 

The genomic DNA of helminths was isolated with a 
OlAamp DNA Micro Kit (Qiagen, Hilden, Germany) 
according to the protocol provided by the manufactur- 
er. DNA concentration and its purity were measured 
using a NanoDrop spectrophotometer (Thermo Scien- 
tific, Waltham, MA, USA). Polymerase chain reactions 
(PCR) were performed on the partial 18S and ITS rDNA 
regions. For nematodes, the ITS region was amplified us- 
ing two sets of primers, Le 18S/26S (Vrain et al., 1992), 
and N93/N94 (Nadler et al. 2005). Amplification of the 
18S region was done with 24F/18P primers (Blaxter 
et al. 1998). DNA of trematodes and cestodes was am- 
plified using 3S/A28 primers for ITS2 region (Bowles 
et al. 1995). PCR were performed in a reaction mixture 
of 20 ul per each sample and consisted of 1 ul of template 
DNA (100 ng of template DNA in 1 ul volume), 1 ul of 
each forward and reverse primer (0.5 uM), 10 ul of 2x 
DreamTaq™ Master Mix (Thermo Fisher Scientific Inc., 
MA, USA) and 7 ul of sterile distilled water. The PCR 
amplifications were performed as described in Nermut’ 
et al. (2015) for the 18S/26S primers, Nadler et al. (2005) 
for N93/N94 primers, Ross et al. (2010b) for 24F/18P, 
and Prasad et al. (2007) for 3S/A28. A volume of 5 ul 
sample of PCR product was run on a 1.5% agarose gel for 
30 min at 100 V to check DNA quality. For rDNA partial 
sequences, PCR products were sequenced with primers 
used for PCR reactions. Sequencing was performed by 
Genomed (Warsaw, Poland). Contigs assembled were de- 
termined using BioEdit version 7.1.3.0 (Hall 1999). 
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Phylogenetic analysis 


The molecular phylogenetic status of helminths was 
determined using BioEdit version 7.1.3.0 for multiple 
sequence alignments (Hall, 1999). Multiple nucleotide 
sequence alignments were generated using also other se- 
quences deposited in GenBank showing the highest sim- 
ilarity to the sequence of examined helminths (Appendix 
I). The species names and GenBank accession numbers 
of the sequences compared to the analyzed helminths are 
shown in the phylograms. Phylogenetic trees were gen- 
erated using MEGA X (Kumar et al. 2018) by maximum 
likelihood (ML) and neighbour-joining (NJ) algorithms. 
The base substitution model was determined for the 18S 
and ITS using MEGA X under the Bayesian Information 
Criterion. The percentage of replicate trees in which the 
associated taxa clustered together in the bootstrap test 
(1000 replicates) are shown next to the branches. Boot- 
strap values above 60% were considered. Branch lengths 
indicate evolutionary rates expressed as the number of 
base differences per site. 


Statistical analyses of prevalence 


A comparison of parasite prevalence was carried out 
between the invasive A. vulgaris and the native A. ater. 
Hybrid slugs were not included since the number of spec- 
imens was low. Generalized linear models (GLMs) with 
binomial distribution were used to test possible differ- 
ences in presence/absence of nematodes and trematodes, 
respectively. All specimens from all investigated popu- 
lations were included in the data set. Locality was used 
as an explanatory factor in addition to slug species to 
test the prevalence of the two different parasite groups. 
statistical analyses were performed in the software R (R 
Core Team 2017). 


RESULTS 
Slug identification 


All collected slugs were identified based on morphology 
and genital morphology and sequences of cytochrome c 
oxidase subunit I (COI) based on BLAST search results. 
In total, 20 European populations were investigated, 
which included 13 populations of A. vulgaris (Norway 
— 9, Poland — 2, France — 1, Germany — 1), and five popu- 
lations of A. ater (all in Norway). Additionally, two pop- 
ulations of hybrids were collected, one population of A. 
ater/A. rufus (Klepp, Norway) and one population of A. 
vulgaris/A. rufus (Zoetermeer, Netherlands). From these 
sites a total of 648 slugs of the genus Arion were col- 
lected comprising 490 A. vulgaris, 92 A. ater, 44 A. vul- 
garis/A. rufus hybrids and 22 A. ater/A. rufus hybrids 
(Table 1). 
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Helminth identification and prevalence 


Helminths were found associated with slugs at 17 of the 
18 sample sites in Europe (94.4%). A total of 501 (77.3%) 
of 648 examined slugs were infected with helminths. All 
slug taxa were infected with nematodes. Trematodes 
were found in A. vulgaris and A. ater, and cestodes in 
A. vulgaris and A. ater/A. rufus hybrids. In A. vulgaris, 
359 specimens were infected with helminths, 1.e., 148 
with nematodes (A//oionema appendiculatum Schneider, 
1859; Angiostoma sp., Phasmarhabditis hermaphrodita 
Schneider, 1859), 198 with trematodes (Brachylaiama 
mesostoma Rudolphi, 1803) and 13 with cestodes (Skrj- 
abinia sp.). In A. ater, 74 of the 92 dissected specimens 
were infected with helminths, 1.e., 57 with nematodes 
(A. appendiculatum, P. hermaphrodita, Angiostoma sp. 
and Entomelas sp.), and 17 with trematodes (B. mesosto- 
ma, Eurytrema sp.). In A. vulgaris/A. rufus, 43 of the 44 
dissected slugs were infected with nematodes (A. ap- 
pendiculatum). In the 22 dissected A. ater/A. rufus, 25 
helminths were found, 1.e., 23 nematodes (4. appendicu- 
latum, P. hermaphrodita), and 2 cestodes (Skrjabinia sp.) 
(Table 1). 

Nematodes were found in slugs at 14 of the 18 visited 
sites (Norway — 11, France — 1, Netherlands — 1, Poland 
— 1; 77.8% of all sample sites), trematodes at 12 of the 
18 sites (Norway — 11, France — 1; 66.7% of all sample 
sites), and cestodes at five of the 18 sites (Norway — 4, 
Germany — 1; 27.8% of all sample sites). A total of four 
nematode species (A. appendiculatum, Angiostoma Sp., 
P. hermaphrodita, Entomelas sp.), two trematode species 
(B. mesostoma, Eurytrema sp.) and one cestode (Skrja- 
binia sp.) species were identified based on morpholog- 
ical and molecular identification (Table 1; explanation 
regarding the three species Entomelas sp., Eurytrema sp., 
and Skrjabinia sp. being listed as “genus sp.” 1s provided 
in the paragraph “Phylogenetic analysis” below). 

From sites positive for nematodes, A. appendiculatum 
was recorded in 11 sites, P. hermaphrodita in eight sites, 
Angiostoma sp. 1n five sites, and Entomelas sp. in only 
one site. For sites positive for trematodes, B. mesostoma 
was found in 11 sites, and Eurytrema sp. in only one site. 
All cestodes were identified as Skrjabinia sp. (Table 1). 


Phylogenetic analysis 


In total, 501 sequences of the 18S and ITS rDNA regions 
were generated (Table 1). These sequences represented 
seven species of helminths from seven families. Sequenc- 
es of the same species were identical across the 18S and 
ITS rDNA regions, so only one representative sequence 
was submitted for each taxon. Obtained sequences of hel- 
minth species have been deposited in GenBank with the 
following accession. numbers: KY355082-K Y355088 
(Table 2). 
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BLAST search for three species sequences did not 
show exact matches with other sequences deposited in 
GenBank. The sequence of K Y355086 revealed 92.98% 
identity with Entomelas dujardini Maupas, 1916 (ac- 
cession number: KF99959]), sequence of KY355087 
revealed 90.46% identity with Eurytrema pancreaticum 
Janson, 1889 (accession number: KY490000), and se- 
quence of KY355088 revealed 83.63% identity with 
Skrjabinia cesticillus Molin, 1858 (accession number: 
AY382321). Other sequences were characterised by 
very low query coverage, resulting in low total scores 
of BLAST searches. None of the sequences available 
in GenBank were significantly similar to sequences ob- 
tained in this study for these three species. Therefore, 
these three detected helminths (1.e., Entomelas sp., Eu- 
rytrema sp., and Skrjabinia sp.) are listed in the study as 
“genus sp.” 

Trees that were inferred from maximum-likelihood 
(ML) and neighbor joining (NJ) revealed identical to- 
pologies. Thus only maximum-likelihood results are pre- 
sented along with bootstrap support from each method 
of analysis (Figs 1—6). The molecular phylogenetic trees, 
generated from partial 18S and ITS of rDNA regions with 
ML and NJ algorithms, showed that the detected nema- 
todes belonged to the species A. appendiculatum (Fig. 1) 
and P. hermaphrodita (Fig. 2), and to the genus Angiosto- 
ma sp., and Entomelas sp. The molecular phylogenetic 
analysis revealed that these two nematodes are closely 
related to Angiostoma margaretae and A. norvegicum 
(Fig. 3), and Entomelas dujardini (Fig. 4). BLAST 
search for Angiostoma sequences revealed 100% iden- 
tity with Angiostoma margaretae Ross, Malan, Ivanova, 
2011 (accession number: HQ115062) and Angiostoma 
norvegicum Ross, Haukeland, Hatteland, Ivanova, 2017 
(accession number: KU712560). Due to the fact that not 
enough material was available for molecular work, the 
28S and COI could not be used for the identification of 
these nematodes. Therefore, the nematodes were listed as 
Angiostoma sp. In the study of Singh et al. (2019), partial 
18S and D2D3 sequences from the same DNA material 
were also found to be almost 100% similar to sequences 
from two different species. The evolutionary history of 
A. appendiculatum was inferred by using the ML meth- 
od based on the General Time Reversible (GTR) model, 
and P. hermaphrodita, Angiostoma sp. and Entomelas sp. 
based on the Tamura 3-parameter model. 

The molecular phylogenetic trees, generated from par- 
tial ITS of rDNA regions with ML and NJ algorithms, 
showed that the detected trematodes belong to the spe- 
cies B. mesostoma and to the genus Eurytrema sp. 
(Fig. 5). The molecular phylogenetic analysis revealed 
that this trematode is closely related to E. pancreaticum. 
The ML method based on the GTR model indicate the 
evolutionary history of B. mesostoma and Eurytrema sp. 
The molecular phylogenetic trees, generated from par- 
tial ITS of rDNA regions with ML and NJ algorithms, 
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Table 1. The prevalence (%) of infection of terrestrial slugs with parasite species collected in Europe 
(N = number of slugs examined). 


Locality 


France: 
Metz 


Germany: 
Gauting 


Netherlands: 
Zoetermeer 


Norway: 
Balestrand 


Hana, Sandnes 


Horten 


Jensvoll, Bodo 


Kjenneveien, 
Fredrikstad 


Kristiansand 


Oslo 


Lindhjem, Larvik 


Manstad, Fredrikstad 


Klepp 


Vesteroya, Sandefjord 


Bergen 


Time, Bryne 


Poland: 
Igolomia 


Rzeszów 
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Parasite species 
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Alloionema ——— KJ851587 


mara Neoalloionema tricaudatum KR617921 


Alioionema sp. KP204849 
gung Neoalloionema sp. KX01/496 
Strongyloides procyonis AB205054 
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| 
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Fig. 1. Unrooted maximum likelihood phylogeny of ITS rDNA regions for A//oionema appendiculatum. The scale bar represents 
0.10 substitutions per nucleotide position. Only bootstrap values above 60% are shown. 
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Fig. 2. Unrooted maximum likelihood phylogeny of ITS rDNA regions for Phasmarhabditis hermaphrodita. The scale bar rep- 
resents 0.10 substitutions per nucleotide position. Only bootstrap values above 60% are shown. 


showed that the detected cestodes belong to the genus Statistical analyses of prevalence 

Skrjabinia sp. (Fig. 6). The molecular phylogenetic anal- 

ysis (ML method based on the Hasegawa-Kishino-Yano A tendency of higher prevalence of trematodes was found 
model) revealed that this cestode is most closely related in A. vulgaris populations compared with A. ater popu- 
to S. cesticillus and Raillietina echinobothrida. lations (Fig. 7; GLM, p=0.0618). On the other hand, a 
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Table 2. The accession numbers of examined slug-parasites with NCBI matches. 


Parasite GenBank no. NCBI match Query | Percentage Source 

family/species coverage identity 

Family: 

Alloionematidae 

Alloionema KY355082 Alloionema appendiculatum; 100 98.95 Nermut' et al. (2015) 

appendiculatum KJ851581 

Family: 

Rhabditidae 

Phasmarhabditis KY355083 Phasmarhabditis hermaphrodita; | 100 100 z 

hermaphrodita FJ516761 

Family: 

Angiostomatidae 

Angiostoma sp. KY355084 Angiostoma margaretae; 100 100 Ross et al. (2011) 
HQ115062 
Angiostoma norvegicum; 100 100 Ross et al. (2017) 
KU712560 

Family: 

Rhabdiasidae 

Entomelas sp 

KY355086 Entomelas dujardini; 100 92.98 Tkach et al. (2014) 

KF999591 

Family: 

Brachylaimidae 

Brachylaima mesostoma 

KY355085 Brachylaima mesostoma;, 100 100 Heneberg er al. (2016) 

KT074964 

Family: 

Dicrocoeliidae 

Eurytrema sp. KY355087 Eurytrema pancreaticum; 99 90.46 Su et al. (2018) 
KY490000 

Family: 

Davaineidae 

Skrjabinia sp. KY355088 Skrjabinia cesticillus; 
AY 382321 81 83.63 - 


somewhat lower prevalence of nematodes 1n A. vulgaris 
than in A. ater (Fig. 8) was revealed, although this differ- 
ence was not significant (GLM, p=0.1569). 


DISCUSSION 


This study presents the occurrence of helminths associ- 
ated with two slug species: A. vulgaris and A. ater, as 
well as A. vulgaris/A. rufus hybrids, and A. ater/A. rufus 
hybrids in some parts of Europe, 1.e., France, Germany, 
Netherlands, Norway and Poland. Previous studies on the 
diversity and distribution of slug parasites have focused 
on the presence of nematodes in slugs (Mengert 1953; 
Gleich et al. 1977; Charwat & Davies 1999; Laznik et al. 
2009; Ross et al. 2010b; Ross et al. 2011; Ivanova et al. 
2013; Ross et al. 2016; Singh et al. 2019). In our study, a 
total of seven species of helminths were found to be asso- 
ciated with these slugs including nematodes, trematodes 
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and one cestode species. The nematodes were identified 
to Alloionema appendiculatum, Angiostoma sp., Phas- 
marhabditis hermaphrodita, Entomelas sp., trematodes 
were identified as Brachylaima mesostoma, Eurytrema 
Sp., and the cestode species was identified as Skrjabinia 
sp. We found a tendency for higher prevalence of trem- 
atodes in A. vulgaris compared with the native A. ater. 
However, we did not find the same pattern in terms of 
prevalence for nematodes. 

Slug-parasitic nematodes were found in all organs of 
the body cavity and also on foot muscles, typical of A. ap- 
pendiculatum. The most intensively studied species of 
slug-parasitic nematode of agricultural and horticultural 
crops 1s P. hermaphrodita. In our survey, P. hermaphro- 
dita was found in eight of the 18 sample sites examined, 
and was found to parasitize both A. ater, A. vulgaris, as 
well as A. ater /A. rufus hybrids. Only juveniles of P. her- 
maphrodita were detected in our study. Within the Rhab- 
ditidae family, Phasmarhabditis is the only genus that 


OZFMK 


Anna Filipiak et al. 


Angiostoma margaretae HQ115062 


Angiostoma norvegicum KU712560 
75:86  Angiostoma dentiferum MK274806 
Phasmarhabditis hermaphrodita FJ516/55 
DO OVL Dhasmarhabditis neopapillosa FJ516754 
| Phasmarhabditis californica KM510210 
sass] Angiostoma dentifera FJ516752 
Phasmarhabditis papillosa KM51021 1 
aw Oscheius chongmingensis EU2/3597 
Oscheius insectivora AFQ83079 
Pellioditis marina AF083021 
Angiostrongylus vasorum AJ920365 
Heterorhabditis indica LN611143 
Steinemema feltiae FJ381667 


me. 
0.020 


Fig. 3. Unrooted maximum likelihood phylogeny of 18S rDNA regions for Angiostoma sp. The scale bar represents 0.020 substi- 
tutions per nucleotide position. Only bootstrap values above 60% are shown. 
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Fig. 4. Unrooted maximum likelihood phylogeny of ITS rDNA regions for Entomelas sp. The scale bar represents 0.10 substitu- 
tions per nucleotide position. Only bootstrap values above 60% are shown. 
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nosy Leucochlondium paradoxum KP9036868 
4 Leucochloridium sp. AY258145 
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Fig. 5. Unrooted maximum likelihood phylogeny of ITS rDNA regions for Brachylaima mesostoma and Eurytrema sp. The scale 
bar represents 0.10 substitutions per nucleotide position. Only bootstrap values above 60% are shown. 
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Fig. 6. Unrooted maximum likelihood phylogeny of ITS rDNA regions for Skrjabinia sp. The scale bar represents 0.20 substitutions 
per nucleotide position. Only bootstrap values above 60% are shown. 
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Fig. 7. Boxplot of trematode prevalence (proportionally) in 
populations of Arion vulgaris and A. ater. 


Is considered to be truly parasitic towards slugs. P. her- 
maphrodita 1s known to be capable of killing many spe- 
cies of slugs from several families (Morand et al. 2004; 
Wilson et al. 2012). The nematode has been formulated 
into an effective biological control agent (Tan & Grewal 
2001; Wilson et al. 2015). Infective dauer juveniles (a 
non-feeding survival stage) of P. hermaphrodita seek out 
the host through a combination of both chemotactic and 
chemokinetic responses towards chemical attractants in 
slug feces and mucus from the foot and mantle (Rae et al. 
2006; Hapca et al. 2007). 

Another nematode known to be associated with ter- 
restrial gastropods is Alloionema appendiculatum. This 
nematode has a broad geographical distribution and has 
been found in areas including Europe, Australia and 
North America (Morand et al. 2004; Laznik et al. 2009; 
Nermut’ et al. 2015). In our study, A. appendiculatum 
was found in 11 of the 18 sample sites examined. The 
nematode was found to parasitize A. ater, A. vulgaris, 
and A. ater /A. rufus and A. vulgaris/A. rufus hybrids. A. 
appendiculatum 1s a common juvenile parasite of many 
terrestrial gastropods. This nematode has both parasitic 
and free-living life stages. During the parasitic cycle, 
third-stage juveniles (J3) enter the slug’s body through its 
foot, where the nematodes moult to the fourth-stage Juve- 
nile (J4) which become encapsulated in the pedal muscu- 
lature. These juveniles then exit the slug and moult into 
free-living immature adults (Cabaret & Morand 1990; 
Laznik et al. 2009). In our study we also identified An- 
giostoma sp., present in five of the 18 sample sites exam- 
ined. The nematode was found to parasitize both A. ater 
and A. vulgaris. With the recent discovery of A. gandav- 
ensis, the total number of described species of the genus 
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Fig. 8. Boxplot of nematode prevalence (proportionally) in 
populations of Arion vulgaris and A. ater. 


Angiostoma is now 19, of which 15 are described from 
molluscan hosts (Singh et al. 2019). International surveys 
reveal that molluscan angiostomatids are present in Eu- 
rope, North America, Africa, South-East Asia and New 
Zealand (Ivanova & Wilson 2009; Ivanova & Spiridonov 
2010; Ross et al. 2010a, b; Ross et al. 2011; Ross et al. 
2017). 

The other nematode identified m our study was En- 
tomelas sp. The nematode was only found in one of the 
18 sample sites examined, and only in the native A. ater. 
Entomelas sp. is classified in the family Rhabdiasidae, 
which includes up to 100 known nematode species par- 
asitic in amphibians and reptiles. All share some mor- 
phological characters but the most remarkable feature 
of rhabdiasids is the regular alternation of parasitic and 
free-living generations (heterogony) in their life cycles 
(Tkach et al. 2014). Entomelas dujardini and E. entome- 
las are commonly associated with Anguis fragilis. Ex- 
perimental infection of the slugs Deroceras reticulatum 
(Agriolimacidae) and Arion subfuscus (Arionidae) with 
infective larvae of E. entomelas and E. dujardini has re- 
vealed that both slug species are classed as paratenic (eu- 
paratenic) hosts for these nematode species (Kuzmin & 
Sharpilo 2000). 

Among the detected parasites were also trematodes, 
Le. B. mesostoma and Eurytrema sp. We found B. me- 
sostoma in 11 of the 18 sample sites examined and it 
was found in both A. ater and A. vulgaris. The genus 
Brachylaima contains 72 species that parasitize mam- 
mals and birds as definitive hosts around the world, ex- 
cept Antarctica. Terrestrial gastropods are involved as 
first and second intermediate hosts. They are important 
from a public health point of view, as they cause diseases 


OZFMK 


Helminths associated with terrestrial slugs in some parts of Europe 21 


in humans like hemorrhagic enteritis, diarrhea, inflam- 
mation of the bile ducts, and anemia (Sirgel et al. 2012; 
Suleman & Khan 2016; Valente et al. 2016). 

The trematode Eurytrema sp. was found in one of the 
18 sample sites, in A. ater only. Species of Eurytrema are 
natural parasites of domestic animals (eg. cattle, goats, 
sheep, pigs, dogs) and wild ruminants (such as buffalos, 
camels, deer) as well as monkeys and humans parasit- 
izing pancreatic ducts and bile ducts. Rarely, terrestrial 
snails of various species (e.g., Bradybaena similaris) are 
intermediate hosts for these parasites (Pinheiro & Am- 
ato 1994). These parasites often cause epithelial hyper- 
plasia, hypertrophy of pancreatic ducts, and periductal 
fibrosis that lead to eurytrematosis (Cai et al. 2012; Man- 
ga-González & Ferreras 2014). 

Terrestrial slugs can also be associated with tape- 
worms. One of the most common tapeworm affecting 
poultry systems 1s Davainea proglottina and the inter- 
mediate hosts are gastropods. Tapeworm segments that 
pass through poultry feces are ingested by snails and 
slugs (of the genera Agriolimax, Arion, Cepaea and Li- 
max) and within three weeks a cysticercoid is produced. 
Adult tapeworms are produced in the infected host 8-15 
days after ingestion of an infected snail or slug (Jordan & 
Pattison 1996). 

In our study, the tapeworm Skrjabinia sp. was found in 
five of the 18 sample sites examined. The tapeworm was 
found in A. vulgaris and A. ater/A. rufus hybrids. Skrj- 
abinia 1s a genus of tapeworms that includes helminth 
parasites of vertebrates, mostly of birds. One of the most 
common parasitic platyhelminths in modern poultry fa- 
cilities throughout the world is S. cesticilTus. It 1s a rela- 
tively harmless species among intestinal cestodes in spite 
of a high prevalence. Sometimes called “broad-headed 
tapeworm”, it infects the small intestine of chicken and 
occasionally other birds, such as guinea fowl and tur- 
key, which are generally in close proximity to backyard 
poultry (Kaufmann 1996; Morishita & Schaul 2007). 
Our study showed that, in some cases, a single slug was 
infected with up to two different species of parasites 
(e. by two species of nematodes or by one species of 
nematode and one species of trematode). Two different 
parasites detected from one host was also sometimes 
observed in a study by Singh et al. (2019). Moreover, 
our study confirmed that endoparasitic helminths appear 
to have a broad host range of slug species (Ross et al. 
2010b; Singh et al. 2019). Alloionema appendiculatum 
was reported from all four different slug species, P. her- 
maphrodita from three, Angiostoma sp. from two, B. me- 
sostoma from two, and Skrjabinia sp. also from two. 

The enemy release hypothesis suggests that species 
become invasive due to a lack of enemies in their in- 
troduced areas (Mitchell & Power 2003; Torchin et al. 
2003). However, Colautti et al. (2004) suggested that 
there are strong, enemy-specific effects on host survival 
and that hosts have developed tailored defenses. Thus, 
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one may expect that it 1s the release from specific ene- 
mies that causes direct changes to survivorship, fecundi- 
ty, biomass, or demographic variables that really matters. 
Alternatively, the loss of enemies against which a host 
is well defended would be of little consequence for the 
host populations (Colautti et al. 2004). We found a con- 
trasting pattern in our study, suggesting that the preva- 
lence of trematodes in invasive slugs is higher than in 
native species. The role of trematodes in slugs should be 
further investigated including comparative studies with 
nematodes. 

The anthropogenic spread of slugs may potentially 
lead to a higher degree of mixing of slug populations 
than would occur when spread naturally. This may also 
include the spread of parasites along with their slug 
hosts, which may have implications for the prevalence 
of trematodes and cestodes that use slugs as intermediate 
hosts and thus may become more common in domestic 
mammals such as cattle and sheep. 
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APPENDIX 


Anna Filipiak et al. 


Species name of helminths, GenBank accession numbers and region of sequences used for phylogenetic analyses. 


Fig. number 


Species name 


Alloionema appendiculatum 
Alloionema appendiculatum 
Neoalloionema tricaudatum 
Alloionema sp. 

Neoalloionema sp. 
Strongyloides procyonis 
Strongyloides fuelleborni 
Strongyloides callosciureus 
Strongyloides robustus 
Rhabditophanes sp. 
Steinernema feltiae 

Angiostoma margaretae 
Angiostoma norvegicum 
Angiostoma gandavensis 
Angiostoma dentiferum 
Phasmarhabditis neopapillosa 
Phasmarhabditis hermaphrodita 
Phasmarhabditis hermaphrodita 
Phasmarhabditis hermaphrodita 
Agfa flexilis 

Caenorhabditis elegans 
Angiostoma margaretae 
Angiostoma sp. 

Angiostoma norvegicum 
Angiostoma dentiferum 
Phasmarhabditis hermaphrodita 
Phasmarhabditis neopapillosa 
Phasmarhabditis californica 
Angiostoma dentifera 
Phasmarhabditis papillosa 
Oscheius chongmingensis 
Oscheius insectivora 

Pellioditis marina 
Angiostrongylus vasorum 
Heterorhabditis indica 
Steinernema feltiae 

Rhabdias bufonis 

Rhabdias engelbrechti 
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GenBank 


accession n? 


KJ851581 
KY355082 
KR817921 
KP204849 
KX017496 
AB205054 
AB272235 
AB272229 
AB272232 
KP204851 
AB243439 
MF192968 
MK214816 
MK214815 
MK214814 
FJ516760 
FJ516761 
KM510202 
KY355083 
MK214813 
FJ589007 
HQ115062 
KY355084 
KU712560 
MK214806 
FJ516755 
FJ516754 
KM510210 
FJ516752 
KM510211 
EU273597 
AF083019 
AF083021 
AJ920365 
LN611143 
FJ381667 
KF999593 
MG428406 


Region 


18S-ITS1-5.8S-ITS2-28S 
18S-ITS1-5.8S-[TS2-288 
18S-ITS1-5.8S-ITS2-288 
18S-ITS1-5.8S-ITS2-288S 
18S-ITS1-5.8S-ITS2-288 
18S-ITS1-5.8S-ITS2-288 
18S-ITS1-5.8S-ITS2-288 
18S-ITS1-5.8S-[TS2-288 
18S-ITS1-5.8S-ITS2-288 
18S-ITS1-5.8S-ITS2-288 
188-ITS1-5.8S-ITS2-285 
18S-ITS1-5.8S-ITS2 
18S-ITS1-5.8S-ITS2 
18S-ITS1-5.8S-ITS2 
18S-ITS1-5.8S-ITS2 
ITS1-5.8S-ITS2 
ITS1-5.8S-ITS2 
ITS1-5.8S-ITS2 
ITS1-5.8S-ITS2 
18S-ITS1-5.8S-ITS2 
18S-ITS1-5.8S-ITS2-288 
18S 

18S 

18S 

18S 

18S 

18S 

18S 

18S 

18S 

18S 

18S 

18S 

18S 

18S 

18S 

ITS1-5.8S-ITS2-28S 
18S-ITS1-5.8S-[TS2-288 


Source 


Nermut et al. (2015) 
This study 

Ivanowa et al. (2016) 
Nermut ef al. (2015) 
Unpublished 

Sato et al. (2006) 
Sato et al. (2007) 
sato er al. (2007) 
Sato et al. (2007) 
Nermut’ et al. (2015) 
Kuwata et al. (2006) 
Unpublished 

Singh er al. (2019) 
Singh er al. (2019) 
Singh er al. (2019) 
Unpublished 
Unpublished 
Tandigan et al. (2014) 
This study 

Singh et al. (2019) 
Imai et al. (2009) 
Ross et al. (2011) 
This study 

Ross et al. (2017) 
Singh er al. (2019) 
Ross et al. (2010b) 
Ross et al. (2010b) 
Tandigan er al. (2014) 
Ross et al. (2010b) 
Tandigan er al. (2014) 
Liu et al. (2012) 
Unpublished (2002) 
Unpublished 

Chilton er al. (2006) 
Unpublished 
Unpublished 

Tkach et al. (2014) 
Kuzmin et al. (2017) 
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Fig. number | Species name — n? Region Source 
Rhabdias bulbicauda KF999600 ITS1-5.8S-ITS2-285 Tkach et al. (2014) 
Rhabdias bermani KF999610 ITS1-5.8S-ITS2-288 Tkach et al. (2014) 
Rhabdias elegans KF999604 ITS1-5.8S-ITS2-288 Tkach et al. (2014) 
Rhabdias pseudosphaerocephala | EU836873 ITS1-5.8S-ITS2-288 Dubey &Shine (2008) 
Rhabdias bakeri EU360831 18S-ITS1-5.8S-ITS2-28S | Dare et al. (2008) 
Rhabdias tarichae MH023523 18S-ITS1-5.8S-ITS2-28S | Johnson et al. (2018) 
Rhabdias picardiae MG195567 18S-ITS1-5.8S-ITS2-28S | Svitin et al. (2018) 
Rhabdias sylvestris KJ018777 ITS1-5.8S-ITS2-288 Tkach et al. (2014) 
Pneumonema sp. KF999603 ITS1-5.8S-ITS2-285 Tkach et al. (2014) 
Pneumonema tiliquae KF99961 1 ITS1-5.8S-ITS2-285 Tkach et al. (2014) 
Entomelas entomelas KF999592 ITS1-5.8S-ITS2-285 Tkach et al. (2014) 
Entomelas sp. KF99960] ITS1-5.8S-ITS2-28S Tkach et al. (2014) 
Entomelas dujardini KF999591 ITS1-5.8S-ITS2-288 Tkach et al. (2014) 
Entomelas sp. KY355086 ITS1-5.8S-ITS2-28S This study 
Steinernema feltiae AB243439 18S-ITS1-5.8S-ITS2-28S | Kuwata et al. (2006) 

5 Leucochloridium paradoxum KP903688 ITS1-5.8S-ITS2-285 Heneberg et al. (2016) 
Leucochloridium sp. AY258145 ITS1-5.8S-ITS2-285 Casey et al. (2003) 
Leucochloridium vogtianum KP903700 ITS1-5.8S-ITS2-288 Heneberg er al. (2016) 
Leucochloridium perturbatum KP903707 5.8S-ITS2-28S Heneberg et al. (2016) 
Brachylaima sp. JX010634 5.8S-ITS2-28S Unpublished 
Brachylaima mesostoma KT074964 5.8S-ITS2-28S8 Heneberg er al. (2016) 
Brachylaima mesostoma KY355085 5.8S-ITS2-28S This study 
Clinostomum album KU708008 18S-ITS1-5.8S-ITS2-28S | Rosser et al. (2017) 
Clinostomum marginatum KU708007 18S-ITS1-5.8S-ITS2-28S | Rosser et al. (2017) 
Macroderoides sp. HQ680850 18S-ITS1-5.8S-ITS2-28S | Tkach & Kinsella (2011) 
Macroderoides texanus EU850398 18S-ITS1-5.8S-ITS2-28S | Tkach et al. (2008) 
Macroderoides flavus HQ680851 18S-ITS1-5.8S-ITS2-28S | Tkach & Kinsella (2011) 
Dicrocoelium hospes EF102026 5.8S-ITS2-28S Maurelli et al. (2007) 
Lutztrema attenuatum KU563718 5.8S-ITS2-28S Unpublished 
Eurytrema sp. KY355087 5.8S-ITS2-288 This study 
Concinnum ten AB521802 5.8S-ITS2-28S Sato et al. (2010) 
Eurytrema pancreaticum KY490000 18S-ITS1-5.8S-ITS2-28S | Suet al. (2018) 
Opisthorchis viverrini MG797539 5.8S-ITS2-288 Sanpool er af. (2018) 

6 Raillietina echinobothrida MH122787 5.8S-ITS2-28S Unpublished 
Skrjabinia cesticillus AY 38232] 5.8S-ITS2-28S Unpublished 
Skrjabinia sp. KY355088 5.8S-ITS2-288 This study 
Parorchites zederi KP893424 18S-ITS1-5.8S-ITS2-28S | Kleinertz et al. (2014) 
Hymenolepis nana LC589873 ITS2 Banzai et al. (2018) 
Mesocestoides litteratus MH936660 ITS1-5.8S-ITS2-288 Unpublished 
Mesocestoides sp. MH936661 ITS1-5.8S-ITS2-288 Unpublished 
Proteocephalus torulosus AY375549 5.8S-ITS2-28S Scholz et al. (2003) 
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Fig. number | Species name 


Raillietina sp. 
Raillietina beveridgei 
Raillietina dromaius 
Raillietina australis 
Raillietina chiltoni 


Hymenolepis diminuta 


Bonn zoological Bulletin 69 (1): 11—26 


Anna Filipiak et al. 


GenBank e 

accession n? Region 
MK201802 5.8S-ITS2 

AY 382318 5.8S-ITS2-28S 
AY 382320 5.8S-ITS2-288 
AY382317 5.8S-ITS2-28S8 
AY 382319 5.8S-ITS2-288 
MF143799 ITS1-5.8S-ITS2 


Source 


Unpublished 
Unpublished 
Unpublished 
Unpublished 
Unpublished 
Unpublished 
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